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Abstract. We explain how the claims of the KPZ scaling theory are con- 
firmed by a recent proof of Borodin and Corwin on the asymptotics of the 
semi-discrete directed polymer. 
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1 Introduction 



The 1986 Kardar-Parisi-Zhang (KPZ) equation [15] is a stochastic partial dif- 
ferential equation modeling surface growth and, more generally, the motion 
of an interface bordering a stable against a metastable phase. The scaling 
theory is an educated guess on the non-universal coefficients in the asymp- 
totics for models in the KPZ universality class. The scaling theory has been 
developed in a landmark contribution by Krug, Meakin, and Halpin-Healy 
[16]. The purpose of our note is to explain how to apply the scaling the- 
ory to the semi-discrete directed polymer. This model has been discussed in 
depth at the 2010 random matrix workshop at the MSRI and, so-to-speak 
as a spin-off, Borodin and Corwin [5] developed the beautiful theory of Mac- 
donald processes, which provides the tools for an asymptotic analysis of the 
semi-discrete directed polymer. As we will establish, the scaling theory is 
consistent with the results in [5], thereby providing a highly non-obvious 
control check. 

To place the issue in focus, let me start with a simple example. Assume 
as given the stationary sequence Xj, j 6 Z, of mean zero random variables 
and let us consider the partial sums 

n 

Sn = Y. X T ( L1 ) 

3=1 

As well studied, it is fairly common that S n / \fn converges to a Gaussian as 
n — )■ oo, i.e. 

lim F(S n < y/Dy/ns) = F G (s) , (1.2) 

n— s-oo 

where Fq is the distribution function of a unit Gaussian random variable. 
Here Fq is the universal object, while the coefficient D > depends on the 
law P and is in this sense model dependent, resp. non-universal. However, 
using stationarity, D is readily guessed as 

oo 

D= E (^oA,) . (1.3) 

j=-oo 

The KPZ class deals with strongly dependent random variables, for which 
partial sums are of size n 1//3 rather than n 1//2 . Fq is to be substituted by the 
GUE Tracy- Widom distribution function, -F GU e, which first appeared in the 



2 



context of the largest eigenvalue of a GUE random matrix [TTJ [22] . -Pgue is 
defined through a Fredholm determinant as 



F G ue(s) = det(l - P s K Ai P s ) . 



(1.4) 



Here 



Ai is the Airy kernel, 



poo 



K Ai (x,y)= d\ Ai(x + X)Ai(y + A) , 



(1.5) 



Jo 



with Ai the Airy function, and P x projects onto the interval [x, oo). For each 
x, P x K Ai P x is a trace class operator in L 2 (K), hence ( }l,4p is well-defined. 
To determine the scale coefficient is less obvious than in the example above, 
but will be explained in due course. Let me stress that the scaling theory is 
crucial for the proper statistical analysis of either physics [27J [28] or computer 
experiments [1J. Without this input, the comparison with theoretical results 
would be considerably less reliable. 

Our paper is divided into two parts. We first explain the scaling theory in 
the context of a specific class of growth models. In the second part the theory 
is applied to the semi-discrete directed polymer model. The convergence to 
the GUE Tracy- Widom distribution is established in [5], of course including 
an expression for the non- universal scale coefficient. Our goal is to explain, 
how this coefficient can be determined independently, not using the import 
from the proof in [5]. 

2 Scaling theory for the single-step model 

In the single step growth model the moving surface is described by the graph 
of the height function h(t) : Z — > Z with the constraint 



hence the name. The random deposition/evaporation events are modeled 
by a Markov jump process constrained to satisfy (12.11) . The allowed local 
moves are then transitions from h(j,t) to h(j,t) ± 2. The dynamics should 
be invariant under a shift in the /i-direction. Hence the rates for the deposi- 
tion/evaporation events are allowed to depend only on the local slopes. It is 
then convenient to switch to height differences 



\h(j + l,t)-h(j,t)\ = l, 



(2.1) 



Vj {t) = h{j + 1, t) - h(j, t) , Vj (t) = ±1 . 



(2.2) 
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A single growth step at the bond + 1) is given by 

T] ->■ T] j ' j+1 , (2.3) 

where r/ J J+1 is the configuration with the slopes at j and j + 1 interchanged. 
If the corresponding rates are denoted by Cjj + i(r)), depending on the local 
neighborhood of + 1), the Markov generator reads 

Lf(v) = £ c j>j+1 ( V ) - f(v)) • (2-4) 

The slope field rj is locally conserved, in the sense that the sum Y^ b j= a v(j^) 
changes only through the fluxes at the two boundaries a, b. To keep things 
concretely, in the following we consider only the wedge initial condition 

h(j,0) = \j\. (2.5) 

The scaling theory is based on the 

Assumption. The spatially ergodic and time stationary measures of the 
slope process r](t) are precisely labeled by the average density 

p = J^2^Ti5> (2 ' 6) 

\j\< a 

with \p\ < 1. 

Our assumption has been formulated more than 30 years ago. Except for 
special cases, it remains open even today, see the book by Liggett [17] for 
more details. The stationary measures from the assumption are denoted by 
p , as a probability measure on {— 1, 1} Z . 
Given n P one defines two natural quantities: 

• the average steady state current 

j(p) = V P (co,i(v)(Vo ~ Vi)) (2.7) 

• the integrated covariance of the conserved slope field 

Mp) = Ww) - /vM 2 ) • (2.8) 
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We first notice that for long times there is a law of large numbers stating 
that 

h{j,t)~tci>(j/t) (2-9) 

for large j, t with a deterministic profile function 0. In fact is the Legendre 
transform of j in the sense that 

0(y) = sup(j/p-jO>)). (2.10) 

|p|<i 

The argument is based on the hydrodynamic limit for nonreversible lattice 
gases [25], which asserts that on the macroscopic scale the density p(x,t) of 
the conserved field satisfies 

-p( x ,t)+—}{p{x,t))=0 (2.11) 

with initial condition 

**.<»={ -i, x<o°: < 2 - i2 > 

The entropy solution to (j23IJ) . fl2TT2l is indeed given by (FZM . fl2~T0|) . 

The average current can be fairly arbitrary, except for the linear behavior 
near p = ±1. On the other hand is convex up with 0(x) = \x\ for |x| > x c . 
At points, where is either linear or cusp-like, the fluctuations may be 
different from the generic case and have to be discussed separately. We set 

A(p) = -j"(p). (2.13) 

Conjecture (KPZ class). Let y be such that is twice differentiable at y 
with 4>"{y) 7^ and set p = 4>\y), \p\ < 1. If A(p) < oo and A(p) ^ 0, then 

lim P(h( [yt\ , t) - t<f>{y) < -(-±\A 2 ) l/ H 1/3 s) = F GVE {s) . (2.14) 

[•J denotes integer part. Since (p"(y) > 0, because of Legendre transform 
AO) < 0. 

On the scale (— |A^4 2 t) 1/ ' 3 the height fluctuations are governed by the 
Tracy- Widom distribution. XA 2 is the model dependent coefficient which, at 
least in principle, can be computed once p p is available. A = for reversible 
slope dynamics, since j = 0. In that case the fluctuations are of scale t 1//4 
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and Gaussian, see [25], Part II, Chapter 3, for a discussion. For nonreversible 
slope dynamics one can still arrange for A = 0. If j is non-zero, there could 
be isolated points at which A vanishes. At a cubic inflection point of j the 
height fluctuations are expected to be of order (t(logt) 1/2 ) 1/2 [211 15]. 

In our set-up the conjecture has been proved by Tracy and Widom [30] 
for the PASEP. In this case the exchange H — to — h occurs with rate p and 
the exchange — h to H — with rate 1 — p, < p < |. Then \i p is a Bernoulli 
measure, hence A(p) = 1 — p 2 , j(p) = \{2p — 1)(1 — p 2 ), and 

-iAA 2 = i(l-2p)(l-p 2 ) 2 . (2.15) 

The profile function is <p(y) = |(1 — 2p)(l + (?//(l — 2p)) 2 ) for \y\ < 1 — 2p and 
0(y) = M f° r M > 1 ~ 2p. For the totally asymmetric case, p — 0, the limit 
(I2.14p has been established before by Johansson [14] . The PushTASEP falls 
also under our scheme with a proof by Borodin and Ferrari [6] . The scaling 
theory is further confirmed for growth models different from single-step, to 
mention the discrete time TASEP [14] , the polynuclear growth model |22j, 
and the KPZ equation [21 125] . 

The theory of Macdonald process [5] has brought a g-deformed version 
of the TASEP in focus. For us here it is a further example for which the 
non-universal constants can be computed. For the rate in (12 .4p we set 

cj,j+i(v) = |(i - %)(i + m+i)g((nj +l (v)) , (2-16) 

where nj is the number of consecutive — slopes to the left of site j. g(0) = 0, 
g(j) > for j > 0, and g increases at most linearly. The q- TASEP is the 
special case where g(j) = 1 — q> , < q < 1, with the TASEP recovered in 
the limit q — > 0. The slope system maps onto the totally asymmetric zero 
range process for length of consecutive gaps between + slope, denoted by Yj, 
j G Z, Yj = 0, 1... . In the stationary measure the Yj's are i.i.d. and 

f Z(a)- 1 , k = 0, , 

¥(Y = k = { )' } k ,_i k ' 2.17 

{ Z(a) l {[[ j=l g{])) a k , fc = l,2,..., 

where 

oo fc 

^(«)=i+E(n^')) _v ( 2 - 18 ) 

k=l j=l 

and a > such that Z(a) < oo. The translation invariant, time station- 
ary measures for the ?7(t)-process with rates (12 . 1 6[) are stationary renewal 
processes on Z with renewal distribution f)2.17p . f )2.18p . 
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The coefficient A can be computed from 

1 - 

ifc, N bg < 6XP [A £ ^ >« = r(A) ' (2 ' 19) 
i=i 

where the average is over the stationary renewal process with parameter a. 
Then p = r'(A) and A = r"(A) at A = 0. The rate function r is implicitly 
determined by 

r(A) = -A - logz(A) , -i-^(A)Z(az(A))e 2A = 1 . (2.20) 

Z(a) 

p, A are computed by successive differentiations. The result is best expressed 
through G(a) = logZ(a). Then 

|(1 + p) = (l + aG')-\ (2.21) 

A = 4(1 + aG')- 3 a(aG'y = -a(l + p)^- . (2.22) 
The average current is given by 

j(p) = -2(c ,i) Q , j = -a(l + p). (2.23) 

One can use (12.23}) together with (I2.2ip to work out A. But no particularly 
illuminating formula results for the combination XA 2 . 

In addition to (I2.14p there is a second scale, which will play no role here, 
but should be mentioned. Instead of the height statistics at the single point 
[yt\ one could consider, for example, the joint distribution of h(ji, t), h(j 2 , t), 
referred to as transverse correlations. The transverse scale tells us at which 
separation |ji — j 2 \ there are nontrivial correlations in the limit t — > oo. The 
KPZ scaling theory asserts that this scale is 

(l\ 2 At) 2/3 . (2.24) 

The factor 1/2 comes from the requirement that the limit joint distribution 
is the two-point distribution of the Airy process. Corresponding predictions 
hold for the multi-point statistics. Also when considering the two-point func- 
tion of the stationary i](t) process, E(?7o(0)^j(t)) — E(r] (0)) , up to a shift 
linear in t, j has to vary on the scale of (I2.24p to have a nontrivial limit as 
t -+ oo. The scale f l2^4j) is confirmed for the PNG [21], TASEP @], and 
PushTASEP [6 J two-point function in case of step initial conditions, for the 
stationary TASEP [10] and stationary KPZ equation [13], and for TASEP 
and PNG [3 IE] in case of flat initial conditions. 



7 



3 The semi-discrete directed polymer model 

Our starting point is a very particular discretization of the stochastic heat 
equation as 

dZj = Zj-idt + Zjdbj . (3-1) 

Here j G Z, t > 0, and {bj(t),j 6 Z} is a collection of independent standard 
Brownian motions. The analogue of the wedge initial condition is 

Z J (0)=5 3<0 . (3.2) 

Hence Zj(t) = for j < and 

dZj = Zj_idt + Zjdbj , j = 1, 2, . . . , 

dZ = Z dbo . (3.3) 

Let us introduce the totally asymmetric random walk, w(t), on Z moving 
with rate 1 to the right. Denoting by E expectation for w(t) with w(Q) = 0, 
one can represent 

Zj(t) = E (exp [ / d6 1u ( s )(s)]<y tt (t) J -)e t . (3.4) 

Zj(t) is the random partition function of the directed polymer w(t), length 
t, endpoints and j, in the random potential dbj(s)/ds. This model was 
first introduced by O'Connell and Yor [20], see also [T5J EES]- I n the zero 
temperature limit one would maximize over the term in the exponential at 
fixed {bj(s)} and fixed end points, see [12] for an early study. The statistics 
of the maximizer is closely related to GUE and Dyson's Brownian motion 

The height corresponds to the random free energy and we set 



hj(t) = log Zj(t), j > , t > . (3.5) 



hj is the solution to 



dhj = e h >- 1 - h idt + db j (3.6) 
and the slope Uj = hj+i — hj is governed by 

duj = (e"^ - e'^-^dt + db j+1 - dbj , j = 1, 2, . . . , (3.7) 

du = e~"° + dbt - db . (3.8) 
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The slope Uj(t) is locally conserved. 

Somewhat unexpectedly, one can still find the stationary and translation 
invariant measures for the interacting diffusions (13. 7p on the lattice Z. They 
are labeled by a parameter r > and are of product form. The single site 
measure is a double exponential of the form 

fi r {dx) = T(r)- l e- e ~ x e- rx dx , r > . (3.9) 

Averages with respect to \i T are denoted by (•),.. The parameters of the 
scaling theory are now easily computed. We find 

p = ( UQ ) r = -ip(r) (3.10) 

with t/> = r'/r the Digamma function on M + . Note that ip' > 0, ip" < 0, 
and p ranges over R. From (13 .7p the random current is —e~ Uj dt — dbj + i and 
hence the average current 

j = -(e-"°) r = -r. (3.11) 

Finally 

A(r) = (u 2 ) r - (u )l = #(r) . (3.12) 

Since the initial conditions force to be convex down, the signs from 
Section 2 are reversed. In particular, the sup in (I2.10p is replaced by the inf 
and 

0(?/) = hlf (-yp-j(-p)), y>Q. (3.13) 



0(0) = 0, 0" < 0, and has a single strictly positive maximum before 
dropping to — oo as y — > oo. Thus t(f)(y/t) reproduces the singular initial 
conditions for (13. 6p as t — > 0. Moriarty and O'Connell [18] prove that 

lim 1^(kJV) = /(k) (3.14) 

iV— >oo iV 



with 



/(«) = inf (ks - V»(s)) ■ (3.15) 



The scaling theory claims that 



lim~V|(t) = 0(y), y>0. (3.16) 
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Hence, using that if)' > and yn = 1, we have 

My) = -/(«) = (s - j/VOO) = mf - yS) , (3.17) 

in agreement with (13.131) . 

The asymptotic analysis of the height fluctuations is due to Borodin and 
Corwin with the result 

Theorem (5.2.12 of [5]). There exists a k* such for < k* < k it holds 

lim ¥(h n (Kn) - ti/(k) < (-2/"(K))- 1 /3 n i/3 s ) = F GUE ( S ) . (3.18) 

n— >oo 

According to (13. lip . A = — j" > 0. Hence in (I2.14p —A is replaced by A 
and the sign to the right of < is +. To see whether with these changes the 
scaling theory is confirmed, we start from 

h [yt} (t) = tcj>{y) + (|Ay4 2 i) 1/3 £ TW (3.19) 

with £tw a GUE Tracy- Widom distributed random variable, hence 

h n {Kn) = nn(f)(K- 1 ) + (§AA 2 /m) 1/3 £ TW . (3.20) 

Now p = —ip(r(p)) is differentiated as 

1 = --0V, = V"(r') 2 + V'r" . (3.21) 

Since —A = j"(p) = —r"(p) and A(r) = ip'(r), one has 

AA 2 = V'V . (3.22) 

Since y = = —r'(p) and y« = 1, we conclude 

AA 2 k = -</>" (3.23) 

and, since / is the Legendre transform of xp, 

\A 2 K = -j j; , (3.24) 

in agreement with (13.181) . 
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4 Conclusion 

The KPZ scaling theory makes a prediction on the non-universal coefficients 
for models in the KPZ class and has been confirmed for PASEP, discrete 
TASEP, and PNG. We add to this list the semi-discrete directed polymer. 
The corresponding stochastic "particle" model is a system of diffusions, Uj(t), 
with nearest neighbor interactions such that the sums ^2jUj(t) are locally 
conserved. This model has a flavor rather distinct from driven lattice gases. 
Still the long time asymptotics in all models is the Tracy- Widom statistics. 
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